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ABSTRACT 


The  Binary  Decision  Diagram  (BDD)  is  a  very  useful  representation  in  the 
design  and  verification  of  switching  functions.  This  is  due  to  to  its  compactness, 
where  size  is  measured  by  the  number  of  nodes.  In  the  implementation  of  logic 
circuits,  connection  of  sub-functions  is  by  means  of  pass  transistors.  The  delay  time 
for  the  interconnections  is  often  larger  than  the  delay  of  the  decision  logic.  For  that 
reason,  crossings  are  often  more  expensive  than  logic.  Planar  Binary  Decision 
Diagrams  are  therefore  desirable  in  implementing  logic  circuits.  This  paper  presents 
a  method  for  finding  a  planar  Ordered  Binary  Decision  Diagram  (OBDD)  for 
threshold  functions.  The  program  that  implements  the  algorithm  is  written  in 
Borland  C-i-i-.  A  special  case  of  Fibonacci  threshold  function  having  up  to  9 
variables  is  analyzed.  It  is  shown  that  Fibonacci  functions  having  up  to  9  variables 
have  planar  OBDD.  With  this  program,  the  characteristics  of  other  threshold 
functions  are  developed. 
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I .  INTRODUCTION 

Boolean  algebra  is  used  widely  in  computer  science  and 
digital  system  design.  Many  problems  in  digital  logic  design 
and  testing,  artificial  intelligence  and  combinatorics  can  be 
expressed  in  a  sequence  of  operations  on  Boolean  functions. 
However,  the  classical  representation  and  manipulation  of 
Boolean  functions  have  many  shortcomings .  A  variety  of  methods 
have  been  developed  for  representing  and  manipulating  Boolean 
Functions.  Those  based  on  classical  representations  such  as 
truth  tables,  Karnaugh  maps,  or  canonical  Siam  of  products  form 
are  quite  impractical — eveary  function  of  n  arguments  has  a 
representation  of  size  2"  [2]  .  A  more  efficient 

representation  is  the  Binary  Decision  Diagram  (BDD) .  It  has 
several  advantages.  Firstly,  most  commonly  encountered 
functions  have  a  reasonable  representation.  For  instance,  all 
symmetrical  functions  are  represented  by  graphs  where  the 
number  of  nodes  grows  at  most  at  a  rate  proportional  to  the 
square  of  the  number  of  arguments  (n^)  .  The  BDD  therefore  has 
a  more  compact  representation.  Secondly,  the  reduced  form  of 
BDD  is  canonical.  (i.e.,  every  function  has  a  unique 
representation) ,  [2] . 

The  binary  tree  can  be  explained  in  the  form  of  a  Shannon 
expansion; 
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If  we  replace  each  node  by  a  2  input  multiplexer,  we  have  the 
following  network: 


X1 


X2 


X3 


X4 


fOOOO  f0010  ^01 00  f0110  f1000  flOlO  filoo  flllO 

foooi  fooii  foioi  torn  fiooi  fioii  fnoi  fiiii 


Figure  3.  Realization  of  logic  function  using  multiplexers. 


A  BDD  is  a  way  to  represent  a  given  function  using  a 
binary  tree.  A  binary  tree  has  a  root  node  which  is  unique  in 
that  it  is  not  a  child  node  of  any  other  node.  The  root  node 
of  a  BDD  represents  the  given  function  f  (x^^ ,X2 ,  .  .  .Xj^)  .  The 
left  child  of  the  root  represents  the  subfunction 
f  (0,X2,X3. .  .x^)  and  the  right  child  f  (1,X2,X3. .  .Xjj)  .  Similarly, 
the  grandchildren,  great-grandchildren  etc.  represent  the 
stibf unctions  associated  with  X2  and  X3  and  so  on,  until  all 
variables  are  assigned.  The  leaf  nodes,  which  have  no  children 
are  assigned  with  constant  0  and  1. 
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AND  OR 

Figure  4.  The  BDD  for  a  AND  and  OR  function. 

Figure  4  shows  examples  of  the  BDD  for  the  AND  and  the 
OR  function.  In  this  case,  the  ordering  of  the  variables 
makes  no  difference  in  the  BDD  structure.  Ref  [2]  showed  an 
interesting  example  of  a  BDD  which  is  highly  dependent  on 
ordering  of  the  variables,  (See  Figure  5)  The  functions  X]_X2 
+  X3X4  +  XgXg  and  Xj^x^  +  X2X5  +  X3Xg  differ  from  each  other 
only  by  a  permutation  of  variables,  yet  the  first  BDD  has  8 
nodes  while  the  second  has  16  nodes.  From  this  example,  one 
can  see  that  a  poor  choice  of  ordering  can  have  very 
undesirable  results. 
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Figure  5.  Example  of  argument  ordering  dependence. 


Besides  the  difference  in  complexity  in  terms  on  number 
of  nodes  /  the  functions  x^x^  +  X2Xg  +  x^Xg  has  crossings  ^  while 
the  other  does  not .  This  shows  that  BDDs  without  crossing  can 
be  found  by  varying  the  order  of  the  variable,  if  they  exists. 

In  the  LSI  implementation  of  a  network,  crossings  are 
expensive  because  they  require  additional  channels  and 
increase  delay.  In  stibmicron  LSI,  the  delay  in  the 
interconnections  are  comparable  to  the  delay  in  the  logic 
elements.  For  a  complex  network,  the  delay  caused  by 
crossings  can  be  unacceptable.  It  is  therefore  understandable 
that  one  of  the  irt^ortant  goals  in  circuit  design  involving 
BDD's  is  to  eliminate,  or  at  least  reduce  crossings  in  the 
network. 

Sasao  &  Butler  [4]  define  a  restricted  planar  BDD: 


5 


Definition  1:  A  function  has  a  restricted  planar  HDD  if 
th®re  exists  a  BDD  without  crossing's,  where  the  1  edge  emerges 
to  the  right  of  the  node  and  0  edge  emerges  to  the  left  of  the 
node;  and  the  constant  1  is  in  the  right  of  the  constant  0. 

This  definition  will  be  used  to  describe  a  planar  BDD  in 
this  text.  Until  recently,  only  HDDs  up  to  5  variables  have 
been  shown  to  have  no  crossing  [4]  .  The  progress  in  this 
field  is  slow  because  manual  enumerations  of  BDD  are  very 
tedious.  To  facilitate  further  research  related  to  BDDs,  the 
process  of  drawing  BDD,  searching  for  those  without  crossings, 
and  producing  their  orderings,  counting  of  their  nodes  etc. 
must  be  computerized. 

The  development  of  such  a  tool  is  one  of  the  main  aims  of 
this  thesis.  The  application  of  this  tool  is  demonstrated  by 
characterizing  a  class  of  functions.  The  following 
definitions  in  accordance  with  [1]  are  used  in  the  subsequent 
discussion: 

Definition  2:  A  switching  function  f(xl,x2,  . . . .xn)  is 
a  mapping  of  f:B”->B,  where  B  =  {0,l}. 

Definition  3:  A  threshold  function  f(xl,x2,  . . . .xn) ,  has 
the  property  that  f  =  1  iff  Wj^Xjj+Wj^.j^x„.3^+. . . . +W;^Xi  ^  T,  where 
T  and  are  integers  and  the  logic  values,  0  and  1  of  x^  are 
viewed  as  integers. 

The  value  of  •  •  •  • for  some  assignment 
of  values  xl,x2,  ....and  xn  is  called  the  weighted  sum.  A 
threshold  function  is  completely  specified  by  a  weight - 
threshold  vector,  (Wn,Wn.i,  . . .  ;T)  .  AND  and  OR  are  special 
examples  of  threshold  functions.  Since  any  switching  function 
can  be  realized  by  a  combination  of  AND,  OR  and  inverter,  any 
switching  function  can  be  realized  by  a  network  of  threshold 
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functions.  In  threshold  logic,  in  general,  each  logic  gate 
represents  more  than  AND, OR  or  NOT  does  in  conventional 
switching  theory.  Accordingly,  the  number  of  gate  that 
realize  a  given  function  is  often  less  than  the  number  of 
gates  required  by  the  AND, OR,  NOT  or  other  conventional  gates. 
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II.  REQUIREMENTS  OF  THE  ANALYSIS  TOOLS 

This  thesis  describes  a  tool  to  analyze  BDD's  for 
threshold  functions .  It  can  be  extended  to  evaluate  other 
functions.  The  tool  is  designed  with  the  following 
considerations  in  mind: 

1 .  Inputs 

The  parameters  describing  the  threshold  functions  such  as 
the  vector  . . . .  ,W2,W3^;  T)  will  be  requested  once  the 
program  is  started.  Several  options  are  provided  to  tailor 
the  analysis  to  certain  requirements.  Since  some  analysis 
requires  the  threshold  to  be  varied,  the  option  to  vary  the 
threshold  from  some  minimal  value  to  a  maximum  value  is  also 
provided.  If  the  aim  is  to  search  for  planar  BDDs,  the 
option  for  producing  various  ordering  of  the  variables  in 
lexicographical  order  is  available.  The  option  for  drawing 
the  BDD  and  printing  in  a  text  file  can  be  activated  to  help 
in  visual  analysis  of  BDD's.  However,  for  certain  BDD 
structures,  when  searching  a  threshold  function  of  many 
variables ,  the  number  of  planar  BDDs  can  be  very  large .  A 
preview  of  which  BDD  is  planar  can  be  made  by  requesting  a 
listing  of  planar  BDDs  first. 

2 .  Processing 

The  program  is  designed  to  evaluate  a  threshold 
function  and  provide  the  result  in  the  form  of  a  truth  table 
for  further  BDD  analysis.  For  other  types  of  functions,  the 
routine  to  convert  the  particular  function  to  a  truth  table 
must  be  written.  The  truth  table  can  be  a  large  table 
having  2”  entries,  where  n  is  the  nxamber  of  variables.  Truth 
tables  are  useful  for  checking  or  interfacing  with  other 
P^O£r^3its .  It  is  also  used  in  the  generation  of  the  initial 
BDD. 

The  program  simplifies  the  initial  BDD  to  a  form  called 
reduced  ordered  BDD.  To  allow  drawing  of  the  BDD,  the  reduced 
BDD  is  laid  out  in  the  form  of  a  binary  tree. 
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At  present,  the  program  is  able  to  analyze  HDDs  of  up  to 
9  variables.  If  threshold  functions  of  more  variables  are 
needed,  the  program  can  be  expanded.  The  program  is  designed 
to  detect  any  crossings  in  the  BDD.  This  feature  would  help 
in  searching  for  planar  HDDs  in  any  type  of  function,  and  is 
useful  for  the  design  of  fast  logic  circuit.  An  important 
output  is  the  nvimber  of  nodes  of  a  reduced  order  BDD.  This 
feature  is  useful  since  the  compactness  of  a  BDD  is  given  by 
the  nxunber  of  nodes  in  the  BDD. 

The  program  produces  all  the  possible  orderings  of 
variables  to  search  for  planar  BDD,  if  it  exists  at  all.  As 
discussed  before,  some  orderings  yield  planar  BDDs. 

3 .  Output 

The  program  provides  the  following  output  for  BDD 
analysis: 

•  Layout  for  drawing  of  the  BDD 

•  Number  of  nodes  for  each  BDD 

•  List  of  BDD's  without  crossings,  showing  the  ordering 
for  each  BDD 
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III.  ALGORITHM  &  IMPLEMENTATION  IN  C++ 

The  algorithm  to  enumerate  binary  decision  diagrams  and 
analyzing  them  is  described  below: 

1.  Generate  Lexicographical  Ordering 

This  module  is  required  for  the  enumeration  of  BDDs.  It 
produces  all  orderings  of  the  variables.  If  two  variables 
have  the  same  weight,  only  the  first  ordering  is  used;  all 
subsequent  orderings  of  variables  corresponding  to  this  same 
arrangement  of  weights  are  ignored.  For  example,  with  weight - 
threshold  vector  (3,2, 1,1;  5),  ordering  X4X3X2Xj^  is  ignored, 
and  only  ordering  X4X3X3^X2  is  used.  The  program  run  time  can 
be  very  long  as  the  number  of  variables  increases  since  there 
are  n!  threshold  functions  to  enumerate. 

2 .  Enumeration  of  Threshold  Function 

An  important  part  of  the  program  is  to  generate  a 
threshold  function  with  weight- threshold  vector 
(  Wj^,  .  .  . ,  W2,Wj^;T  )  which  is  stored  in  an  array  to  be  used  for 
binary  tree  generation  at  a  later  stage. 

3 .  Generation  of  a  Binary  Tree 

A  binary  tree  is  generated  with  2  nodes.  Each  of  the 
nodes  of  the  binary  tree  has  the  following  properties: 

•  Index:  Each  node  is  assigned  an  index  that  increases 
from  top  to  bottom  and  from  left  to  right.  This 
identification  is  used  later  to  check  whether  a  node 
can  be  merged.  Figure  6  shows  an  example  of  the 
indexing  of  the  nodes  for  a  binary  tree  that  describes 
a  4  variables  function.  For  every  node  of  index  i,  its 
left  branch  has  an  index  of  2i,  and  its  right  branch 
has  an  index  of  2i+l. 

•  Pointer  to  left  branch:  Pointer  that  point  to  the  left 
node  immediately  following  that  node  concerned. 

•  Pointer  to  right  branch:  Pointer  that  point  to  the 
right  node  immediately  following  that  node  concerned. 

•  Data:  A  data  value,  DATA,  associated  with  all  nodes 

other  than  leaf  nodes  (those  with  index  from  2”  to 
2(n+i))  initialized  with  the  character  In  the 
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subsequent  steps,  other  values  such  as  0,1  or  an 
alphabet  "a"  ..."j"  will  be  assigned  to  indicate  the 
size  of  the  subfunction  that  has  been  merged.  A  leaf 
node  corresponds  to  some  assignment  of  values  to  all 
variables.  Each  such  node  is  labeled  by  the  value  of 
the  function  for  that  assignment. 


f 


Figure  6.  Indexing  of  Binary  tree. 

4.  Simplification  of  BDD 

Figure  6  is  also  an  example  of  a  full  binary  decision 
diagram  of  a  function.  Such  a  BDD  has  many  redundant  nodes 
and  can  be  simplified.  The  algorithm  to  do  this  is  as  follow: 

a.  If  the  value  of  a  subfunction,  which  is 
stored  in  DATA,  is  a  constant  0  or  l,  terminate 
that  branch.  The  DATA  of  the  nodes  of  the  binary 
tree  that  proceed  after  the  constant  l  or  0  will  be 
filled  with  blanks  to  show  that  the  corresponding 
branches  are  terminated. 
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b.  If  the  current  subfunction  f^^  is  the  same  as 
one  already  generated  fj  ,  move  the  current  branch 
f^  (fj)  over  to  the  other  branch  f  j .  For  example,  if 
fOl  is  the  same  as  flO,  move  flO  branch  over  thus 
merging  the  2  nodes.  The  BDD  resulting  from  the 
above  simplification  is  called  a  Reduced  Order 
Binary  Decision  Diagram  (ROBDD)  [3] .  Letters  "a" 
...  "j"  are  used  to  indicate  the  size  of  the 
subfunction  merged.  If  the  merging  is  done  for  a 
subfunction  with  only  two  children  and  no 
grandchildren,  then  the  merged  subfunction  will  be 
assigned  a  character  "a".  Similarly,  for  nodes  with 
two  children  and  4  grandchildren  will  be  assigned 
the  character  "b"  if  it  is  a  candidate  for  merging, 
etc.  Such  a  scheme  is  implemented  for  nodes  that 
have  up  to  2®  offsprings.  See  the  example  in  table 
1.  The  merging  of  the  subfunctions  is  done  from 
"bottom  up",  which  means  that  comparison  of 
subfunctions  are  done  first  for  subfunctions  with  3 
nodes.  All  subfunctions  with  3  nodes  are  compared 
and  merged  if  found  to  be  identical.  A  merged 
subfunctions  with  three  nodes  will  be  assigned  with 
character  "a".  After  that  the  next  level  of 
subfunctions  with  7  nodes  will  be  compared.  The 
process  goes  on  until  all  subfunctions  are 
evaluated  for  their  potential  to  merge.  The  size 
of  the  largest  subfunction  compared  is  2’^+!. 


Number  of  nodes  in 

subfunction  merged 

3 

7 

15 

31 

63 

127 

255 

Descriptor 

a 

b 

c 

d 

e 

f 

g 

Table  1:  Example  of  descriptor  for  subfunction. 
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c.  If  two  subfunctions  of  one  node  n  are  the 
same,  i.e.  the  node's  left  branch  nQ  and  the  right 
branch  n^  are  the  same,  then  extend  the  branch  to 
n  down  to  the  nQ  and  n^^  eliminating  n.  (i.e.,  do  not 
use  a  multiplexer  ) . 

5 .  Counting  Nodes 

The  niomber  of  nodes  is  a  useful  measure  of  the 
compactness  of  a  BDD.  C++  has  a  recursive  binary  tree 
traversal  function  that  is  very  convenient  for  this  purpose. 
A  valid  node  has  the  characteristics  of  DATA  =  with  two 
children  branches. 

6.  Check  Crossings 

Inorder  traversal  is  used  in  the  checking  of  crossing. 
The  steps  for  inorder  traversal  are  : 

1.  Traverse  left  subtree  using  an  inorder  traversal. 

2.  Visit  the  root  node. 

3.  Traverse  right  subtree  using  an  inorder  traversal. 
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Figure  7 .  Example  of  inorder  traversal .  The  order  of 
traversal  for  the  above  Binary  tree  is  8-4-9-2-10-5-11-1-3-14- 
7-15. 


Using  the  recursive  C++  function  "inorder  traversal", 
the  DATA  elements  of  each  node  of  the  Reduced  Ordered  Binary 
Diagram  are  checked  consecutively  for  a  transition  of  "0" 
to  "1"  .  Not  all  nodes  have  a  0  or  1.  When  a  0  or  l  is 

encountered,  it  is  recorded  to  check  for  0-1  transition.  For 
a  reduced  order  planar  BDD,  there  should  only  be  one  such 
transition  during  the  complete  traversal  of  the  binary  tree, 
i.e..  All  nodes  preceding  the  last  "0"  must  be  "0",  or  any 

symbol  representing  a  merged  subfunction.  Likewise,  any  node 
after  the  first  occurrence  of  a  "1"  must  be  either  "l"  or 
or  any  of  the  symbol  representing  a  merged  siibf unction.  A  BDD 
with  more  than  one  0-1  transition  in  a  inorder  traversal  of 
the  reduced  binary  tree  will  exhibit  crossings,  since  all 
leaves  with  "0"  can  only  be  joined  to  each  other  with 
crossing,  and  likewise  those  leaves  with  "1".  Figures  8 
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Figure  10.  (Step  3.)  Printing  of  simplified  BDD  layout. 


Figure  11.  (Step  4.)  Completion  of  Planar  BDD 

The  current  implementation  of  the  algorithm  is  able  to 
analyze  threshold  functions  of  up  to  9  variables.  Figures  12- 
15.  are  examples  of  6,7,8  and  9  variable  Fibonacci  threshold 
functions  represented  in  the  form  of  planar  BDD. 


Figure  14 .  Planar  BDD  for  f 
vector  (21,13,8,5,3,2,1,1;  34) 


7 .  Output 

There  are  two  main  output  forms: 

a.  Binary  Decision  Diagram:  The  nodes  of  the  planar  BDD 
is  laid  out  in  the  form  of  a  binary  tree  in  ASCII 
text.  For  a  large  tree  (i.e.  of  7  or  more  variable 
threshold  function) ,  the  tree  should  be  viewed  or 
printed  out  in  compressed  font  due  to  the  limited  space. 
The  nodes  can  be  joined  together  in  the  same  way  as  a 
binary  tree.  The  merged  subfunctions  are  represented  by 
alphabets  a,b,....j.  'a'  corresponds  to  the  subfunction 

with  two  leaf  nodes,  b  correspond  to  a  subfunction  with 
two  children  nodes,  each  with  two  leaf  nodes  etc.  These 
nodes  marked  by  a  letter  (a,b, . . . , j)  should  be  joined  to 
the  equivalent  svibfunction  that  it  represent  on  the  left 
of  the  binary  tree.  All  nodes  with  leaf  nodes  of  '0'  and 
should  be  joined  to  sink  nodes  '0'  and  '1'  at  the 
lowest  layer  of  the  binary  tree. 

b.  Listing:  The  listing  of  the  planar  HDDs  without 
crossing  are  printed  to  a  ASCII  text  file  in  the  same 
directory  as  the  program.  Other  characterization,  such 
as  the  number  of  nodes,  ordering  etc.  are  also  listed. 
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IV. 


ANALYSIS 


A.  PLANAR  BDD  FOR  CLASSICAL  THRESHOLD  FUNCTIONS 

Muroga  [6]  defines  canonical  function  as  a  function  of  n 
variables  with  >  W2  >...>Wj^.  He  lists  all  canonical 

positive  threshold  functions  with  five  or  fewer  variables. 

For  six  variables,  only  canonical  positive  self  dual 
functions  are  shown.  A  dual  of  the  function  f  (x)  is  defined 
as 

f'^(x)  =  f'  (X') 

where  f'  is  the  complement  of  f  and  x'  =  (x^^',  X2',...Xjj')  and 
A  function  is  self  dual  if 
f‘^{x)  =  f(x) 

x^'  is  the  complement  of  x„. 

Using  the  methodology  described  in  this  thesis,  the 
ordering  that  produces  planar  HDDs  for  each  of  the  listed 
threshold  functions  is  found.  The  orderings  that  correspond 
to  these  planar  BDD  are  listed  in  Appendix  (A)  .  From  this,  we 
can  state 

Theorem  1  :  All  threshold  functions  of  up  to  five 

variables  and  all  six  variables  canonical  positive  self  dual 
threshold  functions  have  a  planar  BDD. 

Based  on  the  ordering  suggested  by  the  program  written, 
the  BDD  can  be  either  derived  manually  or  by  means  of  the  BDD 
program  whose  source  code  is  listed  in  APPENDIX  (C) ,  for  all 
threshold  functions  listed  in  the  APPENDIX  (A) .  This  table 
greatly  simplifies  the  design  of  efficient  logic  network 
based  on  planar  BDD  for  common  threshold  functions.  Giving 
the  number  of  nodes  that  corresponds  to  each  planar  BDD  helps 
in  further  optimizing  the  logic  design,  by  indicating  which 
ordering  has  the  minimal  nxamber  of  nodes. 
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Theorem  2.  All  threshold  ftmctions  of  up  to  five 
variables  and  all  six  variables  canonical  positive  self  dual 
threshold  fvinction  have  a  planar  BDD  with  the  minimum  number 
of  nodes. 

This  shows  that  planar  BDDs  are  amongst  the  most 
efficient  representations  of  the  threshold  function  with  a 
small  number  of  variables.  This  is  shown  to  be  the  case  for 
threshold  functions  up  to  8  variables. 

The  percentage  of  planar  BDD  decreases  rapidly  as  the 
nvimber  of  arguments  increases,  which  means  that  planar  BDD  are 
extremely  rare  as  the  niomber  of  arguments  increase.  The 
orderings  for  planar  BDD  cannot  be  easily  derived 
heuristically  for  threshold  functions  of  many  variables. 
There  is  often  more  than  one  minimal  BDD.  The  following  table 
shows  the  frequency  of  planar  BDD  versus  all  possible  BDD's; 


No.  of  variables 

1 

2 

3 

4 

5 

6 

Total  no.  of  BDD's 

1 

2 

7 

92 

2578 

1310 

Average  no.  of  BDD's 
that  are  planar 

1 

2 

2.3 

2.7 

4.8 

5.4 

Average  no.  of  BDD's 
that  are  planar  and 

minimal 

1 

2 

1.6 

2.5 

3.86 

4.08 

Percentage  planar 

BDD  {  %  ) 

100 

100 

71.4 

44.5 

16.1 

4.9 

Table  2  :  Planar  BDD  is  very  rare  as  the  number  of  variables 
increases.  Column  6  refers  to  6  variable  canonical  positive 
self  dual  threshold  functions. 
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The  total  number  of  permutations  refers  to  the  total 
number  of  permutations  on  n  variables  that  produce  unique 
functions,  for  all  n-variables  threshold  functions  listed  in 
Appendix  A.  The  average  number  of  permutations  that  are 
planar  is  found  by  dividing  the  total  number  of  planar  BDDs 
by  the  niomber  of  threshold  functions .  When  such  BDD  has 
minimum  number  of  nodes,  it  is  classified  as  planar  and 
minimal . 


Percentage  of  BDD  that  are  Planar 


Figure  17.  Percentage  of  planar  BDD  in  all  unique  permutations. 
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Figure  17  is  a  graph  of  percentage  of  permutations  that  yields 
planar  BDD  versus  the  number  of  variables .  Planar  BDDs  for 
threshold  functions  of  more  than  8  variables  are  very  rare. 
Although  exhaustive  search  planar  BDDs  for  all  threshold 
function  with  8  variables  is  not  done,  search  on  a  small 
sample  of  8  variable  BDD,  such  as  the  one  shown  in  Figure  16 
shows  that  percentage  of  planar  BDDs  approaches  zero. 
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B .  FIBONACCI  FUNCTIONS 


Definition  4:  A  Fibonacci  number  F^  is  specified  by  a 
recursive  relationship;  F^  =  F^.^^  +  •  with  P;|l=P2=^* 

Definition  5;  A  Fibonacci  function  is  a  threshold 
function  with  weight- threshold  vector  . . .  .F2,Tj^}T) , 
where  is  the  i-th  Fibonacci  number  and  0  <  T  <  P„^2* 

A  Fibonacci  function  is  a  threshold  function  with  weight - 

threshold  vector  - F2,Fi;T),  where  F^  is  the  i-th 

Fibonacci  number  and  0  <  T  <  Fn+2-  The  following  have  been 
observed  for  Fibonacci  threshold  functions  of  up  to  9 
variables . 

Theorem  3 .  All  Fibonacci  threshold  functions  with 
thresholds  of  one  to  the  largest  threshold  and  with  an 
ascending  ordering  of  variedsle  have  planar  BDD.s 


This  is  also  stated  in  [l]  .  The  largest  threshold  is 

the  same  as  the  largest  weighted  sum,  1+1+2+3+. _ +Fn=^n+2* 

Table  1  in  Appendix  B  shows  the  number  of  nodes  of  Fibonacci 
functions  for  threshold  of  l  to  the  largest  threshold. 

Theorem  4.  Amongst  the  orderings  that  yield  the  smallest 
niunber  of  nodes  are  (1.)  descending  and  (2.)  ascending  order 
of  weight. 

When  the  threshold  is  1,  the  ordering  does  not  matter  at 
all,  since  this  function  corresponds  to  a  OR  function.  When 
the  threshold  is  at  its  maximum  value,  the  planar  BDD 
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correspond  to  a  flipped  image  of  the  BDD  with  threshold  of  1. 
It  can  be  seen  from  the  listing  in  Appendix  B  that  these 
always  correspond  to  the  minimal  BDD. 

Theorem  5.  The  compactness  profile  of  the  planar 
Fibonacci  HDD's  with  thresholds  from  1  to  its  largest 
threshold  is  symmetrical. 

This  symmetry  can  be  seen  clearly  graphically  as  shown  in 
Appendix  B.  Figure  18  shows  an  example  of  Fibonacci  functions 
of  9  variables,  with  the  plot  of  the  number  of  nodes  versus 
the  threshold  value  of  1  to  their  maximum.  The  distance  from 
one  minimum  to  the  next  one  in  the  compactness  profile 
exhibits  an  interesting  characteristic,  i.e.. 

Distance  of  i^^  trough  to  (i+1)^*^  trough 

- - - — -  ~  Golden  ratio 

Distance  of  (i+l)’-^  trough  to  (i+2)*^^  trough 

Fn 

with  Golden  ratio  =  lim  -  =  1.618 

n-  >  F  ^ 

where  Fj^  is  a  Fibonacci  nvimber  . 
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'll 


Figure  19 .  Distribution  of  Fibonacci  function  by  Nodes  and 
Variables . 

Figure  19  shows  the  distribution  of  nodes  in  the  BDD's  of 
Fibonacci  functions  as  enumerated  by  the  BDD  program.  The 
number  of  variables  and  the  number  of  nodes  in  the  BDD's  are 
plotted  horizontally,  while  the  nximber  of  Fibonacci  functions 
is  plotted  vertically. 


V.  CONCLUSION 


An  effective  method  of  finding  planar  BDDs  for  threshold 
functions  is  developed.  The  algorithm  is  implemented  for 
threshold  functions  of  having  up  to  9  variables.  The  source 
code  in  Borland  C++  is  shown  in  Appendix  C.  It  is 
demonstrated  that  Fibonacci  threshold  functions  having  up  to 
9  variables  have  planar  BDDs .  The  Fibonacci  threshold 
function  is  also  characterized  using  the  algorithm.  It  is 
also  found  that  all  threshold  functions  of  having  up  to  6 
variables  listed  in  [6]  have  planar  BDDs.  The  ordering  that 
produce  planar  BDDs  are  catalogued  in  Appendix  A  for  easy 
reference . 

Follow  up  research  in  this  area  may  produce  an  algorithm 
to  convert  a  BDD  to  its  sum-of -product  expression.  This  would 
facilitate  research  in  more  effective  representation  of 
switching  functions.  Currently,  the  program  does  not  accept 
other  functions  as  input.  It  can  be  made  more  versatile  if 
input  in  the  form  of  a  sum-of-product  expression  to  the  above 
described  program  can  be  implemented. 
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APPENDIX  A.  ORDERING  FOR  PLANAR  BINARY  DECISION  DIAGRAM 

The  following  tables,  extracted  from  [6],  list  all 
threshold  functions  of  3  to  6  variables.  It  is  found  that  all 
the  threshold  function  have  planar  BDD.  Many  of  the  orderings 
are  symmetrical  due  to  identical  weights.  Only  one  of  the 
symmetrical  BDD  is  listed  here.  In  the  following  table,  in 
the  ordering  of  321  means  X3X2X1,  with  X3  as  the  root  node,  and 
the  Xj  as  sink  node  in  the  BDD.  For  each  threshold  function  in 
the  table,  a  disjunctive  form  is  expressed  with  the 
subscripts  of  the  variables;  for  example,  lv23  denotes 
x^.OR.  (X2.AND.X3}  . 

The  planar  BDD  for  the  following  three  variable  threshold 
functions  are  found  and  tabulated  with  their  orderings  and  the 
number  of  nodes,  which  indicates  their  compactness. 


S/N 

Weights; 

Threshold 

Ordering 

Number  of 
nodes 

Threshold  Functions 

1. 

111, -2 

321 

6 

12vl3v23 

2  . 

211, -3 

321 

5 

12vl3 

■■ 

132 

5 

211, -2 

231 

5 

lv23 

■■ 

[imsBBsn 

'"^1  tifST 

132 

r  '  J - -  1 _ ^ 

5 

Table  1.  t’lanar  BDD  tor  three  variable  threshold  function  . 


S/N 

Weights; 

Threshold 

Ordering 

Number  of 
nodes 

Threshold  Functions 

2211;5 

4321 

6 

123vl24 

2431 

6 

1234 

6 

2  . 

2211;2 

4321 

6 

Iv2v34 

2431 

6 

2143 

6 

3. 

1111, -3 

any 

8 

123vl24vl34v234 

order 

4. 

1111; 2 

any 

8 

12vl3v23vl4v34v24 

order 

5. 

2  111;  4 

4321 

7 

123vl24vl34 

1432 

7 

35 


■ 

4321 

1324 

7  . 

2211;4 

4321 

2431 

2134 

8. 

2211;3 

4321 

2431 

2143 

9. 

3211;  5 

4321 

2431 

1432 

1243 

10. 

3211;3 

4321 

2431 

1432 

1243 

11. 

3221;5 

3421 

3241 

1432 

1342 

12  . 

3221;4 

3421 

3241 

1432 

1342 

13  . 

2111; 3 

4321 

1324 

14. 

3111;4 

4321 

1324 

15. 

3111; 3 

4321 

1324 

Table  2 .  Planar  BDD  for  four 


7 

7 


Iv23v24v34 


7 

12vl34v234 

8 

7 

7 

12vl3v23vl4v24 

8 

7 

5 

12vl34 

5 

5 

5 

6 

Iv23v24 

6 

6 

6 

7 

12vl3v234 

7 

7 

7 

7 

12vl3v23vl4 

7 

7 

7 

8 

12vl3vl4v234 

8 

6 

12vl3vl4 

6 

6 

lv234 

6 

variables  threshold  functions 
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S/N 

Weight ; 
Threshold 

Ordering 

#nodes 

Threshold  function 

1. 

22211;2 

54321 

7 

Iv2v3v45 

35421 

7 

32541 

7 

32154 

7 

2. 

mil;  4 

any  order 

10 

1234vl235vl245vl345v2 

345 

3. 

mil;  2 

any  order 

10 

12vl3v23vl4v24v34vl5v 

25v35v45 

4. 

22m;6 

54321 

8 

1234vl235vl245 

15432 

8 

21543 

8 

5. 

22m;2 

54321 

8 

Iv2v34v35v45 

25431 

8 

21543 

8 

6. 

2mi;5 

54321 

9 

1234vl235vl245vl345 

15432 

9 

7  . 

2mi;2 

54321 

9 

1  v2  3  v2  4  v3  4  v2  5  v3  5  v4  5 

15432 

9 

8. 

33211;8 

54321 

7 

123vl245 

35421 

7 

25431 

7 

21543 

7 

9 

33211;3 

54321 

7 

Iv2v34v35 

34521 

7 

25431 

7 

23541 

7 

21543 

7 

21354 

7 

10. 

22211; 6 

54321 

9 

123vl245vl345v2345 

35421 

10 

23541 

10 

32154 

9 

11. 

22211;3 

54321 

9 

12vl3v23vl4v24v34vl5v 

35421 

10 

25v35 

32541 

10 

32154 

9 

37 


32211;7 


32211;3 


33221;8 


33221;4 


mil;  3 


18.  43221;4 


19.  32221;7 


20.  32221;4 


54321 

8 

123vl245vl345 

35421 

9 

32541 

8 

13254 

8 

54321 

8 

Iv23v24v34v25v35 

35421 

9 

32541 

8 

15432 

8 

13452 

9 

13452 

9 

45321 

9 

123vl24vl345v2345 

43521 

9 

24531 

10 

12453 

9 

45321 

9 

12vl3v23vl4v24v34vl5v 

43521 

9 

25 

24531 

10 

21543 

9 

21354 

9 

any  order 

11 

123vl24vl34v234vl25vl 

35v235vl45v245v345 

45321 

8 

123vl24vl345 

43521 

8 

25431 

8 

24531 

8 

14532 

8 

14352 

8 

12543 

8 

12453 

8 

45321 

8 

Iv23v24v34v25 

43521 

8 

25431 

8 

24531 

8 

14532 

8 

12543 

8 

12453 

8 

43521 

9 

123vl24vl34v2345 

43251 

9 

15432 

9 

12543 

9 

43521 

9 

12vl3v23vl4v24v34vl5 

43251 

9 

15432 

9 

14532 

9 

21. 


33111, -7 


54321 

25431 

21543 


7 

7 

7 


33111;3 

54321 

7 

Iv2v345 

25431 

7 

21543 

7 

22111;5 

54321 

10 

123vl24vl25vl345v2345 

25431 

11 

21543 

10 

22111;3 

54321 

10 

12vl3v23vl4v24vl5v25v 

25431 

11 

345 

21543 

10 

33221;7 

45321 

9 

123vl24vl34v234vl25 

43521 

9 

24531 

9 

21543 

9 

21453 

9 

33221;5 

45321 

9 

12vl3v23vl4v24v345 

43521 

9 

24531 

9 

21543 

9 

21453 

9 

33222;7 

54321 

10 

123vl24vl34v234vl25v 

21543 

10 

135v235vl45v245 

33222;6 

54321 

10 

12vl34v234vl35v235v 

21543 

10 

145v245v345 

22211, -5 

54321 

9 

123vl24vl34v234vl25v 

25431 

10 

135v235 

32541 

10 

32154 

9 

22211;4 

54321 

9 

12vl3v23vl45v245v345 

25431 

10 

23541 

10 

32154 

9 

32111;6 

54321 

9 

123vl24vl25vl345 

25431 

9 

15432 

9 

12543 

9 

32111;3 

54321 

9 

Iv23v24v25v345 

25431 

9 

15432 

9 

12543 

9 

39 


33  . 

43221; 8 

43521 

10 

123vl24vl34vl25v2345 

25431 

10 

24531 

10 

14532 

10 

14352 

10 

12543 

10 

34. 

43221;5 

43521 

10 

12vl3v23vl4v24vl5v345 

25431 

10 

24531 

10 

14532 

10 

14352 

10 

12543 

10 

35. 

43321;8 

34521 

10 

123vl24vl25vl35vl34v 

32541 

9 

234 

32451 

9 

15432 

9 

14532 

9 

36. 

43321;6 

34521 

10 

12vl3vl4v23v234v345 

32541 

9 

32451 

9 

15432 

9 

14532 

9 

37. 

43322;8 

32541 

11 

123vl24vl25vl34vl35vl 

15432 

11 

45v234v235 

38. 

43322;7 

32541 

11 

12vl3vl45v234v235v245 

15432 

11 

v345 

39. 

33111;6 

54321 

8 

12vl345v2345 

25431 

10 

21543 

8 

40. 

33111;4 

54321 

8 

12vl3v23vl4v24vl5v25 

25431 

10 

21543 

8 

41. 

22111;4 

54321 

9 

12vl34v234vl35v235vl4 

25431 

12 

5v245 

21543 

9 

42 

33211;6 

54321 

8 

12vl34v234vl35v235 

35421 

8 

25431 

10 

23541 

10 

21543 

8 

21354 

8 

43.  33311;5 


44.  53221;9 


45.  53221, -5 


32211; 6 


32211, -4 


32221;6 


32221;5 


43111;7 


54321 

35421 

25431 

23541 

21543 

21354 


45321 

43521 

25431 

24531 

14532 

14352 

12543 

12453 


45321 

43521 

25431 

24531 

14532 

14352 

12543 

12453 


35421 

32541 

32451 

15432 

13542 

12543 


35421 

32541 

32451 

15432 

13542 

12543 


43521 

43251 

15432 

14532 


43521 

43251 

15432 

14532 


54321 

25431 

12543 


12vl3v23vl45v245 


123vl24vl25vl34 


123vl24vl34vl25vl35v2 

345 


12vl3v23vl4vl5v245v34 

5 


123vl24vl34v234vl25vl 

35vl45 


12vl3vl4v234v235v245v 

345 


12vl345 
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51.  43111;4 


52.  54221; 9 


53.  54221;  6 


54.  54321;9 


55.  54321;7 


54321 

25431 

12543 


43521 

25431 

24531 

43251 

24351 

15432 

14532 

14352 

15243 

12543 


43521 

25431 

24531 

43251 

24351 

15432 

14532 

14352 

15243 

12543 


34521 

25431 

24531 

23541 

23451 

15432 

14532 

13542 

13452 

12543 


34521 

25431 

24531 

23541 

23451 

15432 

14532 

13542 

13452 

12543 


42 


56. 

54322; 9 

25431 

23541 

15432 

13542 

11 

11 

11 

11 

12vl34vl45vl35v235v23 

5 

57  . 

54322; 8 

25431 

11 

12vl3vl45v234v245v235 

23541 

11 

15432 

11 

13542 

11 

58. 

43311;7 

35421 

8 

12vl3v234v235 

32541 

8 

15432 

8 

13542 

8 

59. 

43311;6 

35421 

8 

12vl3v23vl45 

32541 

8 

15432 

8 

13542 

8 

60  . 

42211;7 

54321 

8 

123vl24vl34vl25vl35 

35421 

9 

54231 

8 

32541 

8 

15432 

8 

13542 

9 

61. 

42211; 4 

54321 

8 

Iv23v245v345 

35421 

9 

54231 

8 

32541 

8 

15432 

8 

13542 

9 

62  . 

21111;4 

54321 

11 

123vl24vl34vl25vl35vl 

15432 

11 

45v2345 

63  . 

21111;3 

54321 

11 

12vl3vl4v234vl5v235v2 

15432 

11 

45v345 

64. 

43221;7 

43521 

9 

12vl34v234vl35vl45 

25431 

10 

24531 

10 

24351 

10 

15432 

10 

14532 

10 

14352 

10 

12543 

9 

43 


65. 

43221;6 

43521 

25431 

24531 

24351 

15432 

14532 

14352 

12543 

11 

11 

11 

11 

11 

11 

11 

11 

12vl3vl4v234v235v245 

66. 

32211;5 

35421 

10 

12vl3v234v235vl45 

32451 

10 

15432 

10 

13542 

10 

67. 

31111;  5 

54321 

9 

123vl24vl34vl25vl35vl 

15432 

9 

45 

00 

31111;  3 

54321 

9 

Iv234v235v245v345 

15432 

9 

69. 

53211;3 

54321 

7 

12vl34vl35 

35421 

7 

25431 

7 

23541 

7 

54312 

7 

35412 

7 

15432 

7 

13542 

7 

21543 

7 

12543 

7 

21354 

7 

21345 

7 

70. 

53211; 5 

54321 

7 

Iv23v245 

35421 

7 

25431 

7 

23541 

7 

15432 

7 

13542 

7 

12543 

7 

12354 

7 

71. 

32111;5 

54321 

10 

12vl34vl35vl45v2345 

25431 

11 

14532 

11 

12543 

10 

72  . 

32111;4 

54321 

11 

12vl3vl4v234vl5v235v2 

25431 

11 

45 

14532 

11 

12543 

11 

44 


12vl3v2345 


75.  53321;8 


76.  53321;7 


77.  42111;6 


78.  42111;4 


35421 

25431 

35241 

32541 

23541 

15432 

13542 

15423 

15243 

12543 


35421 

25431 

35241 

32541 

23541 

15432 

13542 

15423 

15243 

12543 


35421 

34521 

32541 

32451 

15432 

14532 

13542 

13452 


35421 

34521 

32541 

32451 

15432 

14532 

13542 

13452 


54321 

25431 

15432 

12543 


54321 

25431 

15432 

12543 


45 


79.  I  42211, -6 


42211;5 


81.  52211;7 


52211;5 


83.  52221;7 


84.  52221;6 


85.  31111;4 


86.  41111; 5 


87.  41111;4 


54321 

35421 

32541 

15432 

13542 

13254 


12vl3vl45v2345 


54321 

9 

35421 

10 

32541 

10 

15432 

10 

13542 

10 

13254 

9 

54321 

7 

35421 

7 

32541 

7 

15432 

7 

13542 

7 

13254 

7 

54321 

7 

35421 

7 

32541 

7 

15432 

7 

13542 

7 

13254 

7 

45321 

9 

43521 

9 

42531 

9 

43251 

9 

15432 

9 

45321 

9 

43521 

9 

42531 

9 

43251 

9 

15432 

9 

54321 

7 

15432 

7 

54321 

7 

15432 

7 

54321 

7 

15432 

7 

12vl3vl4v234vl5v235 


12vl3vl4v2345 


12vl3vl4v234vl5 


12vl3vl4vl5v2345 


12vl3vl4vl5 


lv2345 


Table  3.  Five  variables  threshold  function. 


S/N  Weight; 

Threshold 


332221;7 


222111;5 


3.  433221;8 


4.  332111;6 


5.  322211;6 


6.  543221;9 


433111,-7 


Ordering 


546321 

543621 

216543 

215643 


654321 

365421 

326541 

321654 


326541 

325641 

156432 

154632 


654321 

365421 

265431 

236541 

216543 

213654 


436521 

342651 

165432 

146532 


256431 

254631 

236541 

235641 

156432 

136542 

135642 


365421 

326541 

165432 

136542 

126543 


#nodes 


Threshold  functions 


123vl24vl34vl25vl35vl 

45vl26v234v235v245v34 

56 


123vl24vl34vl25vl35vl 
2  6 vl 3  6 vl 4  5  6 v2  3  4 v2  3  5 v2 
36v2456v3456 


123vl24vl34vl25vl35vl 

45vl26vl36v234v235v24 

56v3456 


12vl34vl35vl36vl456v2 

34v235v236v2456 


123vl24vl34vl25vl35vl 

45vl26vl36vl46v234v23 

56v2456v3456 


12vl34vl35vl45vl36v23 

4v235v2456 


12vl3vl456v234v235v23 

6 
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432211;7 


9.  321111;5 


10.  533211;8 


436521 

265431 

246531 

243651 

165432 

146532 

143652 

126543 


654321 

265431 

165432 

126543 


365421 

245631 

326541 

324651 

146532 

136542 

134652 

165423 


654321 

365421 

326541 

165432 

136542 

132654 


465321 

346521 

432651 

165432 

145632 

143652 


12vl34vl35vl45vl36vl4 

6v234v2356v2456 


12vl34vl35vl45vl36vl4 

6vl56v2345v2346 


12vl3vl45vl46v234v235 

6 


12vl3vl45vl46vl56v234 

5v2346v2356 


12vl3vl4vl56v2345v234 

6 


Table  4.  Planar  BDD  for  6  variables  threshold  functions 


APPENDIX  B.  COMPACTNESS  OF  BDD  FOR  FIBONACCI  FUNCTION 

The  following  table  illustrates  the  compactness  of  the 
BDD  for  fibonacci  threshold  function  of  9, 8,7, 6, 5  variables, 
which  corresponds  to  column  B,C,D,E,F  in  the  table  below.  The 
profile  of  the  number  of  nodes  is  symmetrical  as  the  threshold 
varies  from  1  to  its  maximum.  The  plot  of  theses  nodes  shows 
that  the  ratio  of  one  gap  to  the  next  gap  approximates  the 
golden  ratio.  (Gap  =  distance  from  one  minimal  point  to  the 
next  ) 


Threshold 

1  "■ 

B 

1 

C 

D 

E 

F 

1 

11 

10 

9 

8 

7 

2 

11 

10 

9 

8 

7 

3 

11 

10 

9 

8 

7 

4 

13 

12 

11 

10 

9 

5 

11 

10 

9 

8 

7 

6 

13 

12 

11 

10 

9 

7 

13 

12 

11 

10 

9 

8 

11 

10 

9 

8 

7 

9 

15 

14 

13 

12 

9 

10 

13 

12 

11 

10 

7 

11 

13 

12 

11 

10 

7 

12 

15 

14 

13 

12 

7 

13 

11 

10 

9 

8 

14 

15 

14 

13 

10 

15 

15 

14 

13 

10 

16 

13 

12 

11 

8 

17 

15 

14 

13 

10 

18 

13 

12 

11 

8 

19 

15 

14 

13 

8 

20 

15 

14 

13 

8 

21 

11 

10 

9 

22 

17 

16 

13 

23 

15 

14 

11 

24 

15 

14 

11 

25 

17 

16 

13 

26 

13 

12 

9 

27 

15 

14 

11 

28 

15 

14 

11 

29 

13 

12 

9 

30 

17 

16 

11 

31 

15 

14 

9 

32 

15 

14 

9 

33 

17 

16 

9 

Table  la.  Number  of  nodes  of  fibonacci  function  at 
different  thresholds.  (  To  be  continued  in  table  lb.  ) 
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34  11  10 

35  17  14 

36  17  14 

37  15  12 

38  17  14 

39  15  12 

40  17  14 

41  17  14 

42  13  10 

43  17  14 

44  15  12 

45  15  12 

46  17  14 

47  13  10 

48  17  12 

49  17  12 

50  15  10 

51  17  12 

52  15  10 

53  17  10 

54  17  10 


55  11 

56  17 

57  15 

58  15 

59  17 

60  13 

61  15 

62  15 

63  13 

64  17 

65  15 

66  15 

67  17 

68  11 

69  15 

70  15 

71  13 

72  15 

73  13 

74  15 

75  15 

76  11 

77  15 

78  13 

79  13 


Table  lb  .  (  Cont'd  ) 
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Number  of  nodes 


51 


Number  of  nodes  tq  Number  of  nodes 


APPENDIX  C.  SOURCE  CODE  FOR  BDD  PROGRAM 


/*  Filename:  TREE.H  f 

/*  f 


#define  NULL  0 
typedef  char  DATA; 

struct  node 

{ 

DATA  d; 
int  index; 
struct  node  *left; 

struct  node  *right; 

}; 


typedef  struct  node  NODE; 
typedef  NODE  *BTREE; 


/*********************  Prototypes  ************************/ 

void  print_to_file( int  node_remain,  DATAp[],  BTREE  root,  int 
num_var,  char  filename[],  int  threshold,  int  w[],  int 
x_order [ ] ) ; 

BTREE  ini t_node (DATA  dl,  int  index,  BTREE  pi,  BTREE  p2 ) ; 
BTREE  new_node (void) ; 

BTREE  create_tree (DATA  a[],  int  i,  int  size) ; 
void  inorder (BTREE  root,  int  *cross_flag) ; 
void  preorder (BTREE  root,  int  *node_remain) ; 

void  post  order  (DATA  p[],  BTREE  root,  int  *node_cnt,  int  n\im_var 
)  / 

DATA  merge (BTREE  node_ptr,  int  *node_remain) ; 

void  print_tree (DATA  p[],  BTREE  root,  int  num_var,  int 

x_order [ ] ) ; 

void  cross_test (DATA  d,  int  *cross_flag_ptr) ; 
void  merge_check(  DATA  *f_l[],  int  num_var) ; 

DATA  node (DATA  node_loc [ ] ,  BTREE  root) ; 

void  data_set(int  tf,  int  *threshold,  int  *w_in[]); 

int  sf l_cmp (DATA  sf 1 [ ] ) ; 


55 


/*  Filename:  BDDANATL.C 
/* 

/*  BDD  Tree  Analysis  Program  for  <  10  Variables 
/* 

/*  Functions: 

/*  a.Permutate  all  orderings  of  a  Threshold  Functions 
/*  b.  Evaluate  the  Threshold  Functions 
/*  c.  Generate  BDD  for  each  Function. 

/*  d.  Simplify  the  BDD. 

/*  e.  Examine  for  Crossings. 

/*  d.  Store  results  in  ASCII  file. 

/*  The  results  consist  of 
/* 

/*  i.  Ordering  of  variables  that  produce  planar  BDD 

/*  ii .  Compactness  of  BDD  in  terms  of  no.  of  nodes 

/*  iii.  Binary  Decision  Diagram  of  the  Function 


★  ★  / 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 


#include  "tree.h" 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <math.h> 

#include  <string.h> 

void  main  { ) 

{ 

/*  Counter  for  searching  for  minimum 
no.  of  nodes  */ 

/*  Used  in  evaluation  of  threshold 
function  */ 

/*  Vector  used  for  comparison  of 
sub-functions  */ 

/*  Flag  to  indicate  that  a  valid 
ordering  is  found*/ 

/*  Index  */ 

/*  Index  */ 

/*  Index  */ 

/*  Size  of  subfunction  to  be  compared 


int  min_node=1000 , 

carry [20] , 

cmp_vector [20] , 

found_f lag, 

ii, 

1, 

m, 

sf_size; 


char  sfl[2050],  /*  String  describing  first 

subfunction  to  be  compared  */ 
sf2[2050];  /*  String  describing  second 

subfunction  to  be  compared  */ 


int  t , 


/*  Index  */ 
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k, 
i, 
j  / 

f  [2050] , 
w_in [20], 
w[20]  , 

node_remain=0 , 

*node_cnt , 

nuin_var , 

cross_flag=0, 

x_order [20], 

threshold, 

bin_wt [10]  , 

bin_index[2  050]  , 
tmp_str [2  050]  , 

x[10] [1024] ; 


/*  Index  */ 

/*  Index  */ 

/*  Index  */ 

/*  Enumerated  threshold  function  as 
weighted  sum  */ 

/*  Weights  of  threshold  function 
input  */ 

/*  Weights  of  threshold  functions  */ 
/*  Counter  for  number  of  nodes  in 
binary  tree  */ 

/*  Counter  for  number  of  nodes  left 
in  binary  tree  */ 

/*  Number  of  variable  in  a  threshold 
function  */ 

/*  Indicate  that  there  is  a  crossing 
if  set  */ 

/*  Ordering  of  the  weights  in  a 
threshold  function  */ 

/*  Threshold  in  a  threshold  function 
*/ 

/*  Binary  weights  used  for  computing 
leaf  node  address  */ 

/*  Leaf  node  address  */ 

/*  Temporary  string  used  in  BDD 
manipulation  */ 

/*  Max  number  of  variable  =10  */ 


DATA  f_vector [2050]  , 


input_ok, 
filename [8] , 


/*  Vector  describing  threshold 
function  to  be  converted  into  binary 
tree*/ 

/*  Input  control  flag  */ 

/*  Output  filename  */ 


f_l [2050] , 
f_lo[2050] , 
p[2050]  ; 


/*  Threshold  function  value  in  Is  & 
Os, merging  are  marked  here  */ 

/*  Unchanged  Original  Threshold 
function  vector  */ 

/*  Data  describing  a  node  in  binary 
tree  */ 


FILE  *ofp; 

BTREE  root ; 


/*  Pointer  for  output  file  */ 

/*  Root  address  of  binary  tree  */ 
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/* - Data  Input  - */ 

printf ( "  \n  Number  of  variables  :  " ) ; 
scanf ( "  %d" , &num_var ) ; 

for(i=0;  i<num_var;  i++) 

{ 

printf ( "Weights  :w%d=",i); 
scanf  ("%d",  &;w_in[i]); 

} 

printf ("  \n  Enter  threshold  ="); 
scanf  { "  %d" ,  Scthreshold)  ; 

for(i=0;  i<num_var;  i++) 

{ 

printf ( "Order  X;w%d=  ",i); 
scanf ("%d",  &x_order [ i ] ) ; 

} 

printf ( "Threshold  Function  :  F(  "); 

printf ( "  Weight  :  " ) ; 

for  (  i=0;  i<num_var;  i++) 

printf ("%d  ",w_in[i]); 

printf ( "  ;  %d  ) " , threshold) ; 

printf ( " \n" ) ; 

printf ( "  Order  : " ) ; 
for  (  i=0;  i<num_var;  i++) 
printf ("%d  " , x_order [ i ] ) ; 
printf ( " \n" ) ; 

/*  Loop  for  various  threshold  , provide  the  range  of 
thresholds  that  should  be  evaluated  with  the  febonacci 
functions  here  */ 

for  (threshold=l;  threshold<13 ;  threshold++) 

{ 

ofp=fopen ( " #node . 5ve" ,  "a"); 

f printf (ofp, "Fibonacci  function  with  %d  variables  ,  and 
threshold  =  %d  )  \n" , num_var, threshold  ); 


58 


/******Permutate  for  each  threshold 


for  (i=0;  i<20;  i++) 

{ 

carry [i] =0 ; 
i  f  ( i<nuin_var ) 

x_order  [  i  ]  =nuin_var-i ; 
else 

x_order [i] =0 ; 

} 

x_order [ 0 ] =num_var-l ; 

/★***  generate  increasing  number  for  permutation  generation  */ 
do  { 

x_order [ 0 ] ++ ; 

for  (i=0;  i<num_var;  i++) 

{ 

x_order [ i ] =x_order [ i ] +carry  [ i ] ; 
carry [i] =0 ; 

i f  ( x_order [ i ] >num_var ) 

{ 

x_order [i] =1 ; 
carry [i+1] =1 ; 

} 

}  ■ 

/**start:  ****  filter  out  the  permutation  ******/ 

for  (j=l;  j<=num_var;  j++)  /*  number  to  compare  */ 

{ 

f ound_f lag=0 ; 

for  (i=0;  i<num_var;  i++) 

{ 

if  (x_order [i] ==j ) 

{ 

f ound_f lag=l ; 
break ; 

} 

} 

if  (found_flag==0)  /*  if  found_flag=l  after 

comparing 

all  elements  of  the  array,  it 
is  a  valid  permutation  */ 


break; 

} 


/*********  filter  out  the  permutation  ***********************/ 


if  (found_flag==l) 

{ 

/*  Here  a  permutation  of  an  ordering  is  computated.  The  next 
thing  to  be  done  is  to  find  out  which  of  the  Fibonacci 
function  in  this  ordering  are  planar,  and  has  the  minimum 
number  of  nodes  at  the  same  time  */ 

min_node=1000 ; 

print f ( “ \n\n" ) ; 

printf(''\n  Threshold  =  %d  Ordering  =  ",  threshold  )  ; 
for  (j=0;  j<num_var;  j++) 

printf ( "%d" ,x_order [num_var-j-l] ) ; 
print f ( " \n\n" ) ; 

/***********Pg^j^^^.^|. ends  here  ******************★****/ 


for  (i=0;  i<num_var;  i++) 
w [ i ] =w_in [ (x_order [i ] -1 ) ] ; 


/* - Initialisation - */ 

for  (i=0;  i<10;  i++) 
bin_wt [num_var-i-l] =pow{2 , i) ; 

for  (i=0;  i<(pow(2,10) ) ;  i++) 

bin_index [i] =0 ; 

for  (i=0;  i<num_var;  i++) 

{  for  (j=0;  j<pow  (2 ,  niim_var)  ;  j++) 

x[i]  [ j ] =0; 

} 

for  (i=0;  i<num_var;  i++) 

{  for  (j=0;  j  <pow  ( 2  ,  n\am_var )  ;  j++) 

x[i]  [j  ]  =0; 

} 

for  {j=0;  j<pow(2,num_var+l) ;  j++) 

{ 

f[:)]=0; 

p[j]='  /*  initialise  all  nodes  with  the  symbol 

blanks  */ 

f_l[j]=EOF; 

f_lo[j]=EOF; 

} 
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p[2048] ="\0" ; 

/*  Generate  XO  XI  .  .  .  Xn  where  n=nuin_var - */ 

for  (i=0;  i<nuin_var;  i++) 

{ 

j=0; 
do  { 

for(k=0;  k< (pow(2, i) ) ;  k++  ) 

{ 

x[i]  [j  ]  =0; 
j=j+l; 

} 

for(k=0;  k< (pow (2 , i ) ) ;  k++  ) 

{ 

x[i]  [j]=l; 
j=j+l; 

} 


}  while  (  j<pow(2,nuin_var)  );  /*  while  loop  */ 

}  /*  for  loop  */ 


/* 


generate  results  of  threshold  function  - */ 


for  (  j=0;  j<pow(2,nuin_var)  ;  j++) 
for  (i=0;  i<  nuin_var;  i++) 

{ 

f  [j]=w[i]*x[  (nuin_var-x_order  [i]  )  ]  [j]+f  [j]  ; 


j  ]  =bin_wt  [i]  *x[  (nuin_var-x_order  [i]  )  ]  [  j  ] +bin_index 

j  ]  ; 


[ 


} 


for  (  j  =  0;  j<  (pow(2  ,  nuin_var)  )  ;  j++) 
if  {(f[j]  ==  threshold)  I  I  (f[j]  >  threshold)) 


f_l [ (bin_index [ j ] ) ] = ' 1 ' ; 
f_lo [ (bin_index [ j ])]='!'; 
f  [j]=l; 

} 

else 

{ 

f_l [ (bin_index [ j ] ) ]  = ' 0 '  ; 
f_lo [ (bin_index [ j ] ) ]  = '  0  '  ; 
f  [j]=0; 

} 

} 
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/*  Vector  of  threshold  functions  value  is  formed,  the  next 
step  is  to  simplify  by  comparing  adjacent  subfunctions. 
Identical  subfunction  will  be  combined.  */ 

/* - Adjacent  merge  odd->even - */ 

for  {j=l;  j<num_var;  j++) 

{ 

for  (i=0;  i< (pow (2 , num_var ) ) ;  i++) 

{ 

sfl[i]='\0' ; 

Sf2 [i]  =  '\0'  ; 

} 

sf_size=pow(2, j )  ; 
i  =  0; 

while  ( i< (pow (2 , num_var ) ) ) 

{ 

for  (k=i;  k< (sf_size+i) ;  k++) 

{ 

sfl[k-i]  =  f_lo[k]; 
sf2[k-i]  =  f_lo [ (k+sf_size) ] ; 

} 


if ( (strcmp(sfl,sf2)==0)&&(strcmp(sfl, "OO") !=0) && (strcmp (sfl, 
" 0000" ) 1 =0) && (strcmp (sfl, "00000000" ) ! =0 ) && (strcmp (sf 1 , "00000 
00000000000") !=0) &&( strcmp (sfl, "0000000000000000000000000000 
0000 " ) ! =0) && (strcmp (sfl, "00000000000000000000000000000000000 
00000000000000000000000000000")  !  =0  )&&(  strcmp  (sfl,  "11")  !=0)&Sc 
(strcmp (sfl, "1111" ) ! =0 )&&( strcmp (sfl, "11111111" ) !=0) &&(strcm 
p(sfl, "1111111111111111") !=0) && (strcmp (sfl, "1111111111111111 
1111111111111111" ) !=0) && (strcmp (sfl, "11111111111111111111111 
lllllllllllllllllllllllllllllllllllllllll" )  !-0) ) 

{ 

if ( (sfl !=" 00000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000  ")&&  (sfl!  =  "111111111 111111111111111 11  mil 
111111111111111111111111111111111111111111111111111111111111 
lllllllllllllllllllllllllllllllllllll- ) ) 

{ 

if ( (sfl !=" 00000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
00000000000000000000000000 ")&& (sfl! ="11111111111111111111111 
111111111111111111111111111111111111111111111111111111111111 
111111111111111111111111111111111111111111111111111111111111 
111111111111111111111111111111111111111111111111111111111111 
llllllllllllllllllllllllllllllllllllllllllllllllllllln ) ) 

{ 
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for  (m=  ( i  +  sf_size)  ;  m<  ( i  +  2*sf_size)  ;  in++) 

{ 

if  (sf_size==2) 
f_l  [in]='a'  ; 
else  if  (sf_size==4) 
f_l  [m]  =  'b'  ; 
else  if  (sf_size==8) 
f_l  [m]  = '  c '  ; 
else  if  (sf_size==16) 
f_l [m] = 'd' ; 
else  if  (sf_size==32 ) 
f_l [m] ='e' ; 
else  if  (sf_size==64) 
f_l [m] = ' f ' ; 

else  if  (sf_size==128) 
f_l[m]='g'; 

else  if  (sf_size==256) 
f_l  [rri]  =  'h'  ; 

else  if  ( sf_size==512 ) 
f_l  [rti]  = '  i '  ; 

else  if  (sf_size==1024) 
f_l [m] = ' j  '  ; 

} 

} 

} 

} 

i=i+2*sf_size; 

} 

} 

f_l  [pow  ( 2 ,  niam_var )  ]  = '  \  0  '  ; 

/* - Adjacent  merge  even->odd - */ 

for  (j=l;  j<num_var;  j++) 

{ 

for  (i=0;  i< (pow(2 , num_var) ) ;  i++) 

{  sf 1  [i]  = ' \0 ' ;  /*  Initialisation  */ 

sf2 [i]  =  '\0'  ; 

} 

sf_size=pow (2 , j ) ; 
i=sf_size; 

while  {i< (pow(2 ,num_var) ) ) 

{ 

for  (k=i;  k< ( sf_size+i ) ;  k++) 

{ 

sfl[k-i]  =  f_lo[k]; 
sf2[k-i]  =  f_lo [ (k+sf_size) ] ; 

} 

/*  1=  sf l_cmp (sf 1) ;  */ 
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if ( (strcmp(sfl,sf2)==0)&&(strcmp(sfl, "00") 1 =0) && (strcmp (sfl, 
"0000") 1 =0) && (strcmp (sfl, "00000000") ! =0 ) && ( strcmp ( sf 1 , "00000 
00000000000") l=0)&&(strcmp{sfl, "0000000000000000000000000000 
0000" ) ! =0) && (strcmp (sfl, "00000000000000000000000000000000000 
00000000000000000000000000000" )  !=0)  ScSc  ( strcmp  (sfl,  "11")  !=0)  && 
(strcmp (sfl, "1111") !=0)&&( strcmp (sfl, "11111111") l=0)&&(strcm 
p(sfl, "1111111111111111") !=0)&&( strcmp (sfl, "1111111111111111 
1111111111111111" ) !=0) &&( strcmp (sfl, " lllllllllllllllllllllll 
lllllllllllllllllllllllllllllllllllllllllM )  1=0) ) 

{ 

if ((sfl '="00000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000" ) && (sfl! ="1111111111111111111111111111111 
111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111" ) ) 

{ 

if ( (sfl! ="00000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000 
00000000000000000000000000000000 ")&& (sfl! ="11111111111111111 
111111111111111111111111111111111111111111111111111111111111 
111111111111111111111111111111111111111111111111111111111111 
111111111111111111111111111111111111111111111111111111111111 
11111111111111111111111111111111111111111111111111111111111..) ) 

{ 

for  (m= (i+sf_size) ;  m< ( i+2*sf_size) ;  m++) 

{ 

l=i+sf_size+l ; 

if ( (sf_size==2)  &&  (f_l [1] !  =  'b' )  &&  (f_l  [1] !  =  'c' )  && 

(f_l[l] l-'d')  &&  (f_l[l] !='e')  &&  (f_l[l] !='f')  && 

(f_l [1] !='g' )  &&  (f_l[l] !='h' ) ) 

f_l [m] = ' a ' ; 

else  if  ( (sf_size==4)  && ( f_l [1] !  =  ' c ' )  && ( f_l  [1 ]  !  =  ' d' )  && 

(f_l  [1]  !  =  'e' )  ScSc  (f_l  [1]  !  =  'f  '  )  &&(f  1[1]  !  =  'g' )  && 

(f_l[l] !='h' ) ) 

f_l [m] = 'b' ; 

else  if  ( (sf_size==8)  &&  (f_l [1] !='d' )  &&  (f_l [i] !='e' )  && 

(f_l[l] !='f ')  &&  (f_l [1] !='g' )  &&  (f_l [1] !='h' ) ) 
f_l [m] ='c' ; 

else  if  (  (sf_size==16)  &&  (f_l  [1]  !  =  'e' )  &&  ( f_l  [1]  '  = '  f ' )  && 

(f_l[l] !  =  'g' )  &&  (f_l  [1]  !  =  'h' ) ) 
f_l  [m]  ='d'  ; 

else  if  (  (sf_size==32)  &&  (f_l  [1]  !  =  '  f  ' )  &&  (f_l  [1]  !  =  'g' )  && 

!='h')  &&  (f_l[l] !='i' )  &&  (f_l[l] !=' j ' ) ) 
f_l  [m]  ='e'  ; 

else  if  ( (sf_size==64)  &&  (f_l [1] ! = ' g' )  &&  (f_l [1] ! ='h' )  && 
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(f_l [1]  !  =  'i' )  &&  (f_l [1]  !  =  ' j ' ) ) 
f_l [m] = ' f ' ; 

else  if  ( (sf_size==128)  &&  (f_l [1] !='h' )  &&  ( f_l [1] ! = ' i ' ) 

&& ( f_l [ 1 ] ! = ' j ' ) ) 

f_l[m]='g'; 

else  if ( (sf_size==256) && ( f_l [1] != ' i ' ) & (f_l [1] 1 =' j ' ) ) 
f_l  [m]  =  'h'  ; 

else  if  (  (sf_size==512)&:Sc(f_lo[l]  !  =  ' j  ' )  ) 
f_l  [m]  = '  i  '  ; 

else  if  {sf_size==1024) 
f_l  [m]  =  ' j  '  ; 

} 

} 

} 

} 

i=i+2*sf_size; 

} 

} 

f_l  [pow  ( 2 ,  nuin_var )  ]  = '  \  0  '  ; 
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/*-  Prepare  result  to  be  converted  into  binary  tree  - */ 

for  (i  =  0;  i<  (pow  ( 2  ,  nuin_var ) -1 )  ;  i++) 
f_vector [i] = ' * ' ; 

for  ( i=  (pow (2  ,  nuin_var ) -1 )  ;  i<  (pow (2  ,  (nuin_var+l )  ) +1 )  ;  i++) 

{ 

if  (i==  (pow{2  ,  (nuin_var+l)  )  )  ) 

{ 

f_vector [ i ] = ' \ 0 ' ; 
break ; 

} 

else 

{ 

if  (f_l [i- (pow(2,num_var) -1) ] ==' 0 ' )  /*  Lowest  level 

subfunction  */ 

f_vector [i] = ' 0 ' ; 

if  ( f_l  [i-  (pow(2 , nuin_var)  -1)  ]  == '  1 ' ) 
f_vector [ i ] = ' 1 ' ; 

if  {f_l  [i-  (pow(2,niain_var)  -1)  ]  =='a'  ) 
f_vector [ i ] = ' a ' ; 

if  {f_l  [i  -  (pow(2,niiin_var)  -1)  ]  =='b' ) 
f_vector [ i ] = ' b ' ; 

if  (f_l  [i-  (pow(2 , nuin_var)  -1)  ]  =='c' ) 
f_vect or [ i ] = ' c ' ; 

if  {f_l  [i-  (pow(2,nuin_var)  -1)  ]  =='d' ) 
f_vector [ i ] = ' d ' ; 

if  (f_l  [i  -  (pow(2 , nuin_var)  -1)  ]  =='e' ) 
f_vector [i] ='e' ; 

if  (f_l [i- (pow{2 , num_var) -1) ] =='f ' ) 
f_vector [ i ] = ' f ' ; 

i f  ( f_l  [  i -  ( pow  ( 2 ,  nuin_var )  - 1 )  ]  == '  g ' ) 
f_vector [i] ='g' ; 

if  ( f_l  [i  -  (pow(2 , nuin_var)  -1)  ]  ==  'h' ) 
f_vector [i] = 'h' ; 

if  (f_l  [i-  (pow(2,nuin_var)  -1)  ]  =='i' ) 
f_vector [i] = ' i ' ; 

if  ( f_l  [  i  -  (pow  (2  ,  nuin_var )  -1 )  ]  ==  '  j  ' ) 

/*  Highest  level  subfunction  */ 
f_vector [ i ] = ' j ' ; 

} 

} 

/* - Create  Binary  Tree - */ 

root=create_tree  (f_vector,  0,  (pow(2,  (nuin_var+l)  )-l)  )  ; 

for  (i=0;  i<nuin_var;  i++) 

preorder (root, &node_remain) ;  /*  Call  function  to 

simplify  binary  tree 
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Output 


*/ 


/* 


*/ 


cross_f lag=0 ; 

inorder  (root,  &:cross_flag)  ;  /*  print  out  all  nodes  in 

order  &  do  crossing  check  */ 
node_remain=2 ; 

*node_cnt=2 ; 

/*  Count  the  number  of  nodes  &  free  memory  used  by  binary 
tree--  */ 

postorder (p,  root,  node_cnt,  num_var) ; 


if  {cross_flag==2 ) 

{ 

printf("\n  Planar  BDD\n  " ) ; 
for  (j=0;  j<num_var;  j++) 

fprintf  (ofp,  ■'  %d"  ,x_order  [num_var-j-l]  )  ; 

fprintf (ofp, "  %d  %d\n  ",  threshold,  *node_cnt) ; 
if  (min_node>*node_cnt ) 

J^iri_node=*node_cnt ;  /*  seek  minimxim  nodes  BDD  */ 
^  print_tree (p,  root  ,num_var,  x_order) ; 

if  (cross_flag==4) 

{ 

if  (min_node>*node_cnt ) 

n^i^_riode=*node_cnt ;  /*  seek  minimum  nodes  BDD  */ 
print f (" \nCrossing  occurs  "); 

} 


/*  Write  to  file  if  required - */ 

/*  Ascending  order  */ 
if  (num_var==4) 

{ 

if ( (x_order[0]>x_order [1] ) && (x_order [1] >x_order [2] ) 

&& ( x_order [ 2 ] >x_order [ 3 ] ) ) 

print_to_f ile ( *node_cnt ,  p,  root,  num_var,  filename, 
threshold, w_in,  x_order) ; 

} 

if  (num_var==5) 

{ 
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i  f 

( (x_order [0] >x_order [1] ) && (x_order [1] >x_order [2 ] ) && (x_order [ 
2 ] >x_order [ 3 ] ) && ( x_order [ 3 ] >x_order [ 4 ] ) ) 

print_to_f ile ( *node_cnt ,  p,  root,  num_var,  filename, 
threshold, w_in,  x_order) ; 

} 

if  (num_var==6) 

{ 

i  f 

( (x_order [0] >x_order[l] ) && (x_order [1] >x_order [2] )&&(x_order[ 
2] >x_order [3] ) && (x_order [3 ] >x_order [4] ) && (x_order [4] >x  order 
[5])) 

print_to_f ile ( *node_cnt ,  p,  root,  num_var,  filename, 
threshold, w_in,  x_order) ; 

} 

if  (num_var==7) 

{ 

i  f 

( ( x_order [ 0 ] >x_order [ 1 ] ) && (x_order [1] >x_order [2 ] )&&{x_order[ 
2 ] >x_order [ 3 ] ) && (x_order [3 ] >x_order [4] ) && (x_order [4] >x_order 
[5]  )  &Sc  (x_order  [5]  >x_order  [6]  )  ) 

print_to_f ile ( *node_cnt ,  p,  root,  num_var,  filename, 
threshold, w_in,  x_order) ; 

} 


/*  Descending  order  */ 

if (num_var==4) 

{ 

if ( (x_order [0]<x_order [1] ) && (x_order [1 ] <x_order [2 ] )&&(x_orde 
r [ 2 ] <x_or der [ 3 ] ) && ( x_or der [ 3 ] >x_order [ 4 ] ) ) 

print_to_f ile ( *node_cnt ,  p,  root,  num_var,  filename, 
threshold, w_in,  x_order) ; 

} 

if  (num_var==5) 

{ 

if  (  (x_order  [0]<x_order  [1] )  &&  (x_order  [1  ] <x_order  [2 ]  )Sc&(x_orde 
r [ 2 ] <x_order [ 3 ] ) && {x_order [ 3 ] <x_order [ 4 ] ) ) 

print_to_f ile ( *node_cnt ,  p,  root,  num_var,  filename, 
threshold, w_in,  x_order) ; 

} 

if (num_var==6) 

{ 
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i f ( ( x_order [ 0 ] <x_order [ 1 ] ) && ( x_or der [ 1 ] <x_order [2] ) && { x_orde 
r [2 ] <x_order [3 ] ) && (x_order [3 ] <x_order [4] ) && (x_order [4 ] <x_ord 
er[5])) 

print_to_f ile ( *node_cnt ,  p,  root,  nuin_var,  filename, 
threshold, w_in,  x_order) ; 

} 

if  (n\am_var==7 ) 

{ 

i  f 

( (x_order [0]<x_order[l] ) && (x_order [1] <x_order [2 ] ) && (x_order [ 
2  ]  <x_order  [3  ]  )  &&  (x_order  [ 3  ]  <x_order  [ 4 ]  )  Sc.Sc.  (x_order  [ 4 ]  <x_order 
[ 5 ] ) && (x_order [ 5 ] <x_order [ 6 ] ) ) 

print_to_f ile (*node_cnt ,  p,  root,  num_var,  filename, 
threshold, w_in,  x_order) ; 

} 

}  /*  if  a  permutation  is  found  */ 

}  while  (carry [num_var] ==0) ;  /*  while  permutation 

is  not  completed  */ 

fprintf(ofp,  "\n  Minimal  BDD  has  %d  nodes  \n'' ,  min_node)  ; 
f close (ofp) ; 

}  /*  fibonacci  Loop  ,  from  1  to  T_max  */ 

/*  f close (ofp) ;  */ 

}  /*  end  of  main  */ 


/* - *  / 

/*create  a  linked  binary  tree  from  an  array  */ 

/******************************************** y 

BTREE  create_tree (DATA  a[],  int  i,  int  size) 

{ 

if  (i  >=  size) 

return  NULL; 

else 

return (init_node (a [i] ,  i,  create_tree (a,  2*i+l,  size), 
create_tree(a,  2*i+2,  size) ) ) ; 

} 
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/************************** ! 

/*  Creating  a  binairy  tree  */ 

/************************** y 

BTREE  new_node ( ) 

{ 

return  (malloc (sizeof (NODE) ) ) ; 

} 

BTREE  ini t_node (DATA  dl,  int  i,  BTREE  pi,  BTREE  p2 ) 

BTREE  t; 
t=new_node ( ) ; 
t->index=i+l ; 
t->d=dl; 
t->lef t=pl ; 
t->right=p2 ; 
return  t ; 

} 


/* - Print  out  nodes  in  order  of  left  to  right 

void  inorder (BTREE  root,  int  *cross_f lag_ptr ) 

if  (rooti=NULL) 

{ 

inorder (root->left,  cross_f lag_ptr ) ; 
cross_test (root->d, cross_f lag_ptr ) ; 

/*  print f  ( ''%c“ ,  root  ->d)  ;  */ 
inorder (root->right, cross_flag_ptr) ; 

} 


/* -  scan  for  mergeable  nodes 


void 

{ 


} 


preorder  (BTREE  root,  int  *node_reinain) 

if  (root!=NULL) 

{ 

merge  ( root ,  Scnode_remain )  ; 
preorder (root ->left,  &node_remain) ; 
preorder (root ->right ,  &node_remain) ; 


*/ 


*/ 


/*  count  the  number  of  nodes  left  in  simplified  BDD  - */ 

void  post order (DATA  p[],  BTREE  root,  int  *node  cnt,  int 
num_var)  ~ 

{ 

DATA  rd; 
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int  i,  discount=0; 
if  (root!=NULL) 


postorder  (Sep  [0]  ,  root->left,  node_cnt,  num_var)  ; 

postorder  ( &p  [  0  ]  ,  root->right,  node_cnt,  nuin_var)  ; 

p [root->index] =root->d; 

rd=root->right->d; 

i=root->right->index; 

if  ( (root->d) == ' * ' ) 

*node_cnt  =  *node_cnt+l  ; 

if  (  (rd!  =  '0'  )&&:(rd!  =  'l'  )&&(rd!  =  '*' )  ) 

{ 

if ( (root->left->d=='0' )  I  I (root->left->d== ' 1 ' ) ) 
dis count =1 ; 

if  (  (i>(pow(2,  (nuxn_var-l)  ) -1 )  &&  ( i<pow(2 ,  nurri_var)  )&:&(rd=='a' )  )  ) 
dis count =1 ; 

if  I  (i>(pow(2,  (num_var-2)  ) -1)  &&  (i<pow(2,nuin_var-l)  )&Sc(rd=='b' 
dis count =1 ; 

^  ( 2 ,  (nuin_var-3  ))-!)&&( i<pow  ( 2 ,  nuin_var-2  )  )  &&  (rd== '  c ' 

dis count =1 ; 

if  I  ( i>  (pow  (2 ,  (num_var-4 )  )-!)&&(  i<pow  ( 2 ,  nuin_var-3  )  )  &&  { rd==  '  d' 
dis count =1 ; 

if  I  (i>(pow(2,  (nuin_var-5)  ) -1 )  &&  (i<pow(2  ,  nuin_var-4 )  )&&(rd=='e' 
discount=l ; 

ff)  ^  '  (nunL-var-6 )  )-!)&&(  i<pow  (2 ,  nuin_var-5 )  )  &&  ( rd==  '  f  ' 

discount=l ; 

if  I  (i>(pow(2,  (num_var-7)  )  -1 )  &&  (i<pow(2 ,  nuin_var-6)  )&&(rd=='g' 
dis count =1 ; 

^  '  (num_var-8 )  )-!)&&;(  i<pow  (2  ,  nuin_var-7 )  )  &&  (rd==  '  h ' 

discount=l ; 

if  (  ( i>  (pow  (2 ,  (num_var-9 )  )-!)&&( i<pow (2 ,  nuin_var-8 )  )  &&  ( rd== '  i  ' 
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discount=l ; 

if  (  ( i>  (pow  (2  ,  (nuin_var-10 )  )-!)&&( i<pow  (2 ,  num_var-9 )  )  &&  (rd==  '  j 
')  )  ) 

discount=l ; 
if (discount==l ) 

{ 

*node_cnt=*node_cnt-l ; 
dis count =0 ; 

} 

} 

free (root ) ; 

} 


/**************  simplification  routines 


for  Binary  Tree  *****/ 


/*  Return  a  '0/1'  if  left  branch  &  right  branch  are  both 
'0/1'  or  else  mark  node  with  '*'  to  indicate  that  merging  is 
not  possible  */ 


DATA  merge (BTREE  node_ptr,  int  *node_remain_ptr ) 

{ 

DATA  t , In , rn ; 
ln=node_ptr->lef t->d; 
rn=node_ptr->right->d; 

i  f  ( ( ln==rn ) && ( node_ptr->index 1=1)) 

{ 

if 

( (ln==' 0' )  I  I  (ln=='l' )  I  I {ln=='a' )  I  I  (ln=='b' )  I  I  (ln=='c' )  I  I  (ln= 
=  'd' )  ) 

{ 

t=node_ptr->lef t->d; 
node_ptr->d=t ; 

node__ptr->lef t->d= '  ';  /*  terminate  branch  */ 

node_ptr->right->d= '  '; 
return  t ; 

} 

else  if  ((ln=='e')  M  (ln=='f')  ||  (ln=='g')  ||  {ln=='h') 

I  I  (ln== ' j  ' )  ) 

{ 

t=node_ptr->left->d; 
node_ptr->d=t ; 

node_ptr->left->d='  ';  /*  terminate  branch  */ 
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node  ptr->riqht->d= '  ' 

return  t ; 

} 


} 

} 

/★ -  Print  result  - * 

void  print_tree  (DATA  p[],  BTREE  root,  int  nuin_var, 
x_order [ ] ) 

{ 

int  i , j ; 

print f ( " \n\n" ) ; 

printf ( "X%d 

%c " , x_order [ 0 ] , p [ 1 ] ) ; 

printf ( " \n\n" ) ; 

printf ( "X%d  %c" ,x_order [1] ,p [2] ) ; 

printf ("  %c",p[3]); 

printf ( " \n\n" ) ; 

printf  (■'X%d  %c"  ,x_order  [2]  ,p[4]  )  ; 

for  (j=5;  j<8;  j++) 

printf {"  %c",p[j]); 

if  {nuin_var>2) 

{ 

printf ( " \n\n" ) ; 

if  (nuin_var !  =3  ) 
printf {"X%d  " ,x_order [3] ) ; 

printf { "%c" ,p [8] ) ; 
for  (j=9;  j<16;  j++) 
printf ("  %c",  p[j]); 


if  {nuin_var>3) 

{ 

printf ( " \n\n" ) ; 

if  (nuin_var !=4) 
printf ("X%d  " ,x_order [4] ) ; 

printf ( "%c" ,p[16] ) ; 
for  (j=17;  j<32;  j++) 
printf  (■■  %c",p[j]); 

if  (n\iin_var>4 ) 

{ 


int 
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printf("\n\n  " ) ; 
for  (j=32;  j<64;  j++) 
print f ( "  %c" ,p [ j ] ) ; 

} 

} 

} 

} 

/* - Print  result  to  file - */ 

void  print_to_file{int  node_cnt,  DATA  p[],  BTREE  root,  int 
num_var , 

char  filename [],  int  threshold,  int  w_in[],  int  x_order[]) 

int  i,j; 

FILE  *ofp; 

ofp  =  fopen ( "D_perm. 7vd" , "a" ) ; 

fprintf (ofp, " XnThreshold  function  F= ( " ) ; 
for  (  i=0;  i<num_var;  i++) 
fprintf  (ofp,  nd  '■,w_in[i]); 
fprintf (ofp, " ;  %d  )  \n",threshold); 

fprintf (ofp, "Order  :  " )  ; 
for  (  i  =  0;  i<n\am_var;  i++) 
fprintf (ofp, "X%d  " ,x_order [i] ) ; 

fprintf (ofp, " \nNumber  of  nodes  =  %d  \n\n" , node_cnt ) ; 


fprintf (ofp. 


******  Literal  ***************************** 
"\Binary  Tree  Map\n"); 


/ 


for  (i-64;  i< (pow(2 , (num_var+l ) ) )  ;  i++) 


if  (p[i] 1='  ' ) 

^  fprintf (  ofp,  "Index  =  %d  %c  \n",i,  p[i]); 
fprintf (ofp,  "\n\n"); 


★  ★  ★  ★ 


/ 


/  level  1  */ 

fprintf (ofp, " \n\n" ) ; 
fprintf (ofp, " 

%c",  p[l])  ; 


level  2  */ 
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f print f (ofp, " \n\n" ) ; 
fprintf (ofp, " 
fprintf (ofp, " 

%c",p[3]  )  ; 


if  (nuin_var>l) 

{ 

fprintf  (ofp,  "  \n\n'' )  ; 
fprintf (ofp, " 
for  (j=5;  j<8;  j++) 
fprintf (ofp, " 

} 

if  (nuin_var==2 ) 
fprintf (ofp, " \n\n" ) ; 


%c"  ,p[2]  )  ; 


-level  3  */ 


%c",  p[4]); 


%c",  p[j]  )  ; 


if  (nuin_var>2) 

{ 

fprintf (ofp, " \n\n" ) ; 

/*  if  (nuin_var !  =4 ) 
fprintf (ofp, "X%d 
fprintf (ofp, “  %c  ",  p[8]); 

for  (j=9;  j<16;  j+4-) 
fprintf (ofp, "  %c", 

} 

if  (nuin_var==3 ) 
fprintf (ofp, " \n\n" ) ; 


-level  4  */ 


" , x_order [ 4  ]  )  ; 


c " ,  P  [  j  ]  )  ; 


-level  5  */ 


if  (nuin_var>3) 

{ 

fprintf (ofp, "\n\n" ) ; 
fprintf (ofp, "  %c",  p[16]); 

for  (j=17;  j<32;  j-)-+) 
fprintf (ofp, "  %c",  p[j]); 


if  (nuin_var==4 ) 
fprintf (ofp, "\n\n" ) ; 


if  (nuin_var>4) 


-level  6  */ 
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{ 


fprintf (ofp, " \n\n" ) ; 
fprintf (ofp, "  %c",p[32]); 

for  (j=33;j<64;  j++) 
fprintf  (ofp,  "  %c'’,p[j]); 

} 

if  (nuin_var==5 ) 
fprintf (ofp, " \n\n" ) ; 

/*  level  7  */ 

if  (nuin_var>5) 

{ 

fprintf (ofp, "\n\n" ) ; 
for  (j=64;  j<128;  j++) 
fprintf  (ofp,  "%c  ■■,p[j]); 

fprintf (ofp, "\n\n\n\n" ) ; 

} 

f close (ofp) ; 


/*  Examing  if  there  is  crossing - */ 

cross _ flag  =0  :  no  '0'  or  '1'  occur  yet 

cross_flag  =1  :  first  '0'  occured 

cross_flag  =2  :  first  '1'  occured 

cross_flag  >=4  :  >1  '0->l'  transition  has  occured 

void  cross_test (DATA  d,  int  *cross_flag) 


if  (  (*cross_flag==0)  Sc&  (d==' * ' )  ) 

*cross_f lag=0 ; 

if  ( (*cross_flag==0)&&(d=='0' )  ) 

*cross_f lag=l ; 

if  (  ( *cross_f lag==l )  &Sc  (d== '  1 '  )  )  /*  First  crossing  occurs 

all 

the  time  */ 

*cross_f lag=2  ; 

if  (  (*cross_flag==2)&:&(d=='0' )  ) 

*cross_f lag=3 ; 
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/*  >  1  crossing  occurs 


if  ( (*cross_flag==3) && (d=='l' ) ) 
*/ 

*cross_f lag=4 ; 

} 


/*****■>:*****■)!************************************************ 

Compare  sub-function  to  see  that  sfl  is  not  all  '0'  or  all  '1' 
before  EOF;  return  an  integer  0  if  the  above  is  true,  else 
return  an  integer  1  return  -1  is  error  calling  of  function 

*******************************■),*****************************/ 

int  sf l_cmp (DATA  sfl  [  ] ) 

{ 

int  num_char=0 , j  =  0 ; 
int  cnt01=0; 

while  (sfl [j+1] !=EOF) 

{ 

if  (  (sfl  [j]  ==sfl  [j+1]  )&:&(sfl[j  +  l]  !=EOF)&&(  (sfl  [j]=='0' )  I  |  (sfl 
[j]=='l'  ) ) ) 
cnt01++; 

D++; 

} 

while  (sf 1 [num_char] ! =EOF) 
num_char++ ; 

return  (num_char-cnt01 ) ; 

}  /*  end  of  function  */ 


/* - end - END - END - END - END - */ 
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